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Abstract. We give a hypergeometric proof involving a family of 2-variable Krawtchouk 
polynomials that were obtained earlier by Hoare and Rahman [SIGMA 4 (2008), 089, 
18 pages] as a limit of the 9 ~ j symbols of quantum angular momentum theory, and shown 
to be eigenfunctions of the transition probability kernel corresponding to a "poker dice" 
type probability model. The proof in this paper derives and makes use of the necessary 
and sufficient conditions of orthogonality in establishing orthogonality as well as indicating 
their geometrical significance. We also derive a 5-term recurrence relation satisfied by these 
polynomials. 
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1 Introduction 

It was in the SIDES meeting in St-Adele near Montreal that one of us (MR) presented a paper 
reporting the discovery, by his co-author Michael Hoare and himself [2008], of a "new" system 
of 2-variable Krawtchouk polynomials, orthogonal with respect to a trinomial distribution. The 
motivation of their paper was to find eigenvalues and eigenfunctions of the transition probability 
kernel: 

min(Ji,ji) min(i2,i2) 

KAUi,J2;ii,i2) = ^ ^ b{ki,ii;ai)b{k2,i2;a2) 

ki=0 k2=0 

X hiji -ki,J2 -k2\N -ki -k2;l3i,f32), (1.1) 



where b{x,N;p) = ()p^{l — p)^ ^ is the binomial distribution, while 

62(x, y; N-p,q)=(^ V-(/J'(l -p- qf-^-y, 

is the trinomial, both normalized to 1. In (jl.ip the parameters {ai,a2, f3i,l32) are probabilities 
of a two-step cumulative Bernoulli process, and hence necessarily in (0,1), while {ii,i2) and 
{ji-,32) represent the initial and final states of the process. Borrowing a result from the angular 
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momentum theory of quantum mechanics the authors of [12] were able to show that the 2- 
dimensional Krawtchouk polynomials 

0<i+j+k+l<N 

do indeed satisfy the requirements for them being the eigenfunctions of (|l.ip . where {x,y) 
represents the state-variable and {m, n) the spectral parameters. It goes without saying that 
for (ll.2p to be an orthogonal system with respect to a distribution of the form b2{x, y; N; r]i, 772) 
the parameters ui, vi, U2, V2 must be related to r]i, r/2, as well as satisfy some additional 
conditions among themselves. It was found in [12], again with a cue from the physics literature 
that these conditions are all satisfied provided the u's and v's are parametrized in the following 
way 

(Pi + P2){pi + Ps) {Pi + P2){P2 + Pi) 

Ul = 7 ; ; ; r, U2 — 



PiiPi +P2+P3+ Pa) P2{Pi +P2+P3+ Pa) 

(pi + P3){P4 + P3) {P2 + Pa){P3 + Pa) f. „s 

P3[Pi + P2 + P3 + Pa) Pa[Pi + P2 + P3 + Pa) 

and consequently, 

^ PlP2{Pl +P2+P3+Pa) ^ P3Pa{PI +P2+P3+Pa) ^ ^s 

{Pl + P2){Pl + P3){P2 + PaY ^ {P2 + Pa){P3 + Pa){P1 + P3)' 

For the origin of the work in [12j the reader may consult [3l [9l [TOl [TT] . 

Fortunately, in the audience, a very attentive listener, Masatoshi Noumi, was there to point 
out to (MR) that these polynomials are not only not new, but a special case of the multivariable 
generalization of the Gaussian hypergeometric function: 



n n 



U{-x^)^ Ui-mi)^ 

1=1 j=i 



F;»>(-x,-m;-iV;u) = ^ ^^— •^^^^, (1.5) 






where the a^j's are nonnegative integers taking values from to n, such that ^ctij < N, 



«j 



N being assumed as a nonnegative integer. Here we are following the notation in [T^. In 
the original definition of Aomoto and Gelfand [H [6], A^ need not be an integer, nor the x's 
and m's. Furthermore, the space on which their functions are introduced was a bit more general, 
a generalization we do not need for our purposes. To be sure, these authors' primary interest 
was not to look at (jl.Sp as a multidimensional extension of the Krawtchouk polynomials: 

Pnix)=2Fi{-x,-n;-N;7]~^), (1.6) 



rather some structures that they contain. H. Mizukawa [17j proved that the functions in (jl.Sp 
are the zonal spherical functions on a very special class of Gelfand pairs made up of the reflection 
groups G{r, 1, n) and the symmetric group Sn- For a very nice account of a way to obtain many 
discrete orthogonal polynomials in terms of certain {n + 1,771 + 1) hypergeometric functions see 
also the work of Mizukawa and Tanaka [19]. As mentioned above, we learned from Professor 
M. Noumi that these functions give the multivariable Krawtchouk polynomials independently 
obtained in [12]. In a very recent paper by Mizukawa, see [18], he has established the ortho- 
gonality of Krawtchouk polynomials in n variables by using very different techniques from the 
ones in this paper. See also the additional comment at the end of our paper. 
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The origin of the work in [3l El [T0| \TT\ [T2] is the analysis of a very concrete probabilistic 
model, namely "poker dice". Its corresponding eigenfunctions are seen in [T2] to be given in 
terms of a family of polynomials that, as indicated above, are now identified with the Gelfand- 
Aomoto polynomials. It is likely that this may be the first probabilistic application of the 
Gelfand-Aomoto polynomials. They may also be applicable to other models in the physical 
sciences. 

One should point out that the hypergeometric functions involve both parameters as well as 
variables. Depending on the issue at hand one can consider these functions as depending on 
one or the other set of "variables" . This is already the case in the classical one variable case as 
indicated by the expression (jl.6p above. The fact that these Krawtchouk polynomials (or the 
higher level Hahn polynomials) could be so useful in analyzing naturally appearing models in 
statistical mechanics was not anticipated in the classical book by W. Feller [5], where one can 
read about the Ehrenfest as well as the Bernoulli-Laplace models. For several applications of 
the Krawtchouk polynomials to several parts of mathematics see [16]. For a very good general 
guide to the field see [l]. 

Our objective in this paper is less ambitious in one sense and more in another - namely, that 
we still restrict ourselves to the n = 2 case, but not necessarily on the reflection group, but 
to the general situations where the parameters Ujj's will be determined by the requirement of 
orthogonality. Hoare and Rahman [12] have done that problem, but we will approach it from 
a different angle. We will refrain from parametrizing the Uij right at the outset, instead looking 
for conditions they must satisfy among them in order that the 2-variable polynomials: 

(2) 

F^ '{-mi,-m2]-xi,-X2;-N;ui,vi,U2,V2) 



-E 



{-mi)i+j(-m2)k+i{-xi)i+k{-X2)j+i i j ki 



i\jlkW.{-N)i+,+k+i 
become orthogonal with respect to the trinomial: 



U\v{ulv2 = Pm.i,m2{xi,X2) (1.7) 



\Xl,X2j 

It may be worth mentioning that a prior knowledge of this weight function is not essential 
since one could easily derive it by using the binomial generating function of the polynomials 

-^ m,nl,X, y )■ 

We now state the main results in the paper, namely (|1.8|) and (|1.9p below. 
In Sections [21 [3] and Hj we will show that the necessary and sufficient conditions of orthogo- 
nality are: 

(a) r]iui + r]2Vi = 1, 

(b) riiU2+mv2 = l, (1-8) 

(c) rjiUiU2 + ri2ViV2 = 1, 

with the r/'s assumed to be given such that < ryi, r/2 < 1 and r/i + r/2 < 1. One can easily verify 
that these three conditions are all satisfied by (|1.3p and (II. 4p . 

However, one of the main reasons for going back to this problem is to find a 5-term recurrence 
relation for (II. 7p , which interestingly, is more easily found by using the p's as in (II. 3|) and (II. 4p 
than using (jl.Sp instead. In Section [5] we'll show that, if we denote (jl.7p by Pmi,m2{xi-,X2) then 

[N - mi - 1712) < 7 : TT r {Pmi+l,m2{Xl,X2) - -Pmi.ma (^l , X2)) 

I [Pl + P3){PlP4 - P2P3) 

P2Pi{Pl + P3)(Pl +P2+P?,+Pi) ,jj , s_p , ^^ 

{P2+P4){PlPi-P2Pz) 
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Pi +P3 

-7712 \ {Pmi,m.2-l{xi,X2) - Pmi ,m2 {x , y)) 

P2 +PA 

= {{Pl +P2)xi - (P3 + PA)x2)Pm^,m2{xi,X2). (1.9) 

This recursion relation is vahd when the variables xi, X2 are nonnegative integers taking 
values whose sum is at most A^. 

If we insist on a difference operator in the variables mi, 7712 with an eigenvalue that is linear 
in xi, X2 and involves only the four nearest neighbours of 777i, 77x2 this is essentially the only 
choice. This was the result of extensive symbolic computations carried out beginning at the 
time that [8] was written. This statement was proved in general in [14] a paper that kindly 
acknowledges this work carried out initially in a special case. If one allows the eight nearest 
neighbours we get another linearly independent difference operator, a fact also found by symbolic 
computation by us and independently proved in |14j . 

There is by now a rather large literature dealing with orthogonal polynomials in several 
variables. A reference that is still useful is [13| Vol. 2]. A comprehensive treatment is found 
in [3] . For some of the issues that we are interested in the reader can see \lb\ [7] and the references 
in these papers. 

It may be useful to point out that when the conditions (jl.Sp are not met the polynomials Pm,n 
still satisfy difference equations in the indices {m,n), of the type given in [151 [7j. When the 
conditions (jl.Sp are met these recursions become much simpler in that they involve a smaller 
number of neighbouring indices. Having a recursion relation involving the smallest number of 
neighbours of the index {m, n) might be important in certain numerical implementations of these 
recursions as well as in potential signal processing applications of these polynomials. In those 
situations, having a minimal number of sampling points could be a useful feature. 

In dealing with the same polynomials Iliev and Terwilliger, see the very recent nice paper |14] . 
have found two 7-term recurrence relations. In fact a referee has kindly pointed out that our 
5-term relation ()1.9p can be derived by using a combination of these. 

It is clear that both in |14] and [H] the replacement of the variables Uij by the consideration 
of a matrix with entries 1 — Uij, properly augmented, played a very important role. This matrix 
is considered in [M] and an interpretation is given in terms of Lie algebras. In [18] the author 
proves that the orthogonality of the columns of this matrix with respect to a weight built out 
of the r]i is the appropriate extension of (II. Sp . 

The F^ notation used in the 2-variable case (ll.7p and more generally in the 77-variable 
case (jl.Sp is a reflection of the fact that these are generalizations of the standard Appell-Kampe 
de Feriet function 

F,(a;M';c;...)=EE%f^-V- 
A very useful integral representation of this Fi function is the double integral 

r(c) 



/ / ^i-'^2~\^-^i-^2r"'~''~Hi-^ix-^2yr''d^id^2, 

Jo Jo 



T{b)r{b')r{c-b-b')J, 

provided 

0<Re{b,b',c-b-b'). 

(2) 
This extends to F^ as well, which can be easily verified: 

■p / \ 

Ff'{ai,a2:,bi,b2;c;ui,vi,U2,V2) 



r(ai)r(a2)r(c - ai - 02) 
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Jo Jo 
x{l-viCi-V2a)~'"d^id^2, (1.10) 

which we shall find very useful in our calculations, even though the parameters in the case 
of (jl.7|) do not satisfy the convergence conditions of the integral in (jl.lOp . We will take the 
point of view that whatever identities we find by using (I1.10|) with < Re(ai, 02, c — oi — 02), 
are also valid where oi, 02, c are, in fact, negative integers. 

For the direct hypergeometric proof that we are planning to give in the following pages it 
will be necessary to make use of the transformation formulas: 

(2) 
Fl '{ai,a2;bi,b2;c;ui,vi,U2,V2) (1.11) 

.^ ^_aw, N-a2 77(2)/^ , , , Ui-Vi -Vi U2-V2 -V2 \ 

= (l-t;i) Hl-^'2) ""Fl '[ai,a2;bi,c-bi-b2;c;- ,- ,- ,- 

V 1-^1 ^-Vl l-f2 1-'W2/ 

= {l-ui) '{I-U2) ^Fl '[ai,a2;c-bi-b2,b2;c;- ,- , :; ,- , 

V 1-ui l-ui I-U2 I-U2 J 



which were proved in Hoare and Rahman [12] . But there is a third transformation that we shall 
find occasions to use, that is valid when (mi, 1712) and (rEi,rE2) are pairs of nonnegative integers, 
as is A^, satisfying the triangle inequality: < mi + m2 <N,0<xi + X2<N, and that is 

77(2)^ M \ (Xi +X2- N)rm+m2 /, , „n 

F{ '{-mi,-m2;-xi,-X2;-N;ui,vi,U2,V2) = 7— -^ (1.12) 

[ -'Vjmi+m2 
(2) 

X F{ ' {-mi, -1712; -xi, -X2;N + I - xi - X2 - mi - 7712; 1 - mi,1 - vi,l - U2, 1 - V2), 
which is just a generalization of the transformation: 

2Fi{-m,-x;-N;u) = — — , ^ 2Fi{-m,-x;N + I — x -m;l — u) 

= ^7377r^2i^i(-m, -x;N+l-x-m;l-u). (1.13) 

In fact ()1.12p and (|1.13p easily extend to the multidimensional case F^"" , provided one is dealing 
with terminating series. It may be remarked here that for (J1.12p and (J1.13p to be true, indeed 

(n) 

in the general case of F^ , the parameter A^ need not even be an integer. 



2 A general expression for orthogonality sum and proof of (11.81) 

Let us denote 

CV,m2 = X^X] ^2(2:1, 3;2;iV;r/i,772)-Pmi,m2(3;i,2;2)-Pni,n2(a;i, 2:2). 



Xi,X2 

0<xi+X2<N 

At {nil, 1^2) = (0, 0), and {rii, ^2) ^ (0, 0), this simply represents a generating function for these 
polynomials, namely: 

4T' = (1 - ^1^1 - V2Vir'{i - mu2 - mv2r. 

So, at the three points (0,0), (1,0) and (0, 1) the pairwise orthogonality between the first and 
the last two simply amounts to the conditions (a) and (b) given in (jl.Sp . To obtain the condition 
at the points (1,0) and (0, 1) we need some more computations. 
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For ease of computation we will imagine, for the time being, that —mi, —m2, —N are complex 
numbers ai, 02, c such that < Re(ai, 02, c — ai — 02). Then using p.lO|) we get 



Tic) 



Tni,n2 _ 

'"^''"^ r(ai)r(a2)r(c-ai-a2) 



/ / rfeirf6er"'e2^"'(i-ei-6 

Jo Jo 



ic— ai— a2 — 1 



: Y^ Y^ b2{xi,X2] N; r?i, 7?2)-Pni,n2 (a;i, 3;2)(1 - Ciui - C2U2T' (1 - 6^^i - 6^'2)'''- 



Xl X2 



Clearly 



^^62(2:1, 3^2; A^;??l,??2)(-Xl)j+fc(-2;2)j+/(l -6""! - ^2U2T^il -6^1 - ^2V2T^ 



Xl,X2 



(-7v),+,-+fc+,r?i+s^2^'(i - 6^1 - C2U2y+\i - 6^1 - ^2V2y+' 

X {1 - 6(f?l^il + f?2t'l) - 6(m^2 + ??2W2)}^"*"^"''"' 

(-iv),+,-+fc+K^i(i - emi - e2n2))'+'(r/2(i - 6^1 - C2V2)y+', 



(2.1) 



by virtue of (jl.Sp fa) and (jl.Sp fb). For general {1711,1712) and {111,172), we use ()2.ip and recast 
back to the original parameters, getting 



7-711,712 
-'mi, 7712 



E 



-"-l)i+j(-"'2)fc+« 



i!j!A;!/! 
^ {-i-j-k-l)mi+m2 p(2) 



(??iui)*(r72t^iy(m^2)*^(??2^'2)' (2.2) 

F| '{-mi,-m2;-i-k,-j-l;-i-j-k-l;ui,vi,U2,V2). 



Let us take (7711,7712) = (0,1) and (771,772) = (1,0), so that 

4,f = E ^47i^(mni)^(7?2t^i)^^^^^Fi(-l; ~i, -j; -7 - j;7Z2,7;2) 



^,J 



iljl 



-N) 



(1 '''\Y^^{mui)\mviy{i + j)2F, 



N 



«j 



7!j! 



l-V2sr^ (-l)j-f j(??l^il)'(%'"l)■' 



A^ 



E 



«j 



iljl 



J2F1 



-j)2Fl 


"-l,-7_ 'ii2 
-I -J 1- 


-7^2" 
-7^2 . 


-1,-7 1 - 7i2 






i 


'1-772 







1 '''.Y^^{lnulyil^2v^y[J+^ 



N 



«j 



iljl 



l-V2f , 1-7X2 

7727^1 + 7?i7ii 



N 



1 - 7;2 



1 - 7i2 
1-772, 

(7/i7il(l - 7i2) + ??2t'l(l - V2))/N 



which must vanish, so (jl.Sp fc) must be satisfied in addition to (jl.Sp fa) and (jl.Sp fb). 

It is worth noting that by solving the first two conditions of (|1.8p one can show that the 
third condition amounts to 



U1V2 = U2V1, 



(2.3) 



where 



Ui = l-ui\ 



Vi = l- V- 



1,2. 



Condition (j2.3p has a simple geometrical interpretation as a cone embedded in four dimen- 
sional space. In a subsequent paper we will look at a geometrical interpretation for the corre- 
sponding orthogonality conditions in the case of more than two variables. 
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We would also like to point out that if % and % are parameters of 62(^,1, ^,2; -^;??i,%) for 
the dual orthogonality of the P's, then they must satisfy 

(a) r/iui +772^2 = 1, 

(b) fj^vi+fJ2V2 = l, (2.4) 

(c) ri^uiVi+ri2U2V2 = l. 

3 Reduction of /"^'"^ 

By the first transformation in (|l.lip F^ inside the sum in (j2.2p becomes a multiple of Fi, 
which, transformed by (J1.12p gives 

i,j,k,l 

^^ {-J - Omi+m2 jp f ■ , • , ; , 1 1 - ni 1 - n2 

Fi I -t - k;-mi,-m2;j + I + 1 -mi -1712] 



(--/V)mi+m+2 V 1 - Ul 1 - Z;2 

Set i + k = r,j + l = s,i = r — k, j = s — I, to get 

IP f ,1 ^2-^2 

X i^i -n2; -r, -s; rei + 1 - r - s; — , — 
V ui Vi 

' (-^)™i+™2 jp^ (_^. _^^^ _„^^. , + 1 _ ,„^ _ ^2; 1^^, i^^^l . (3.1) 



(-A')mi+m2 V l-IJl l-t;2 

Since 

IP f ,1 ^2 t'2 

(-ni - n2)r+s „ f ^ U2 ^ V2 

ti -n2; -r, -s; -m - 712; 1 , 1 



(-ni)r+5 V ui vi 

by (fTT^ . (IXTD reduces to 

Tm,n2 (1 - vi)"'^ (1 - ^2)"^^ v^ (-ni - n2)r+. , .^. ., 

-^mV,m2 = TT^Tn 2^ IM (^1^1) (^21^1) 

/ U2 U2 

X Fi -n2] -r, -s; -ni - n2; 1 , 1 

V Ul Vl 

( \ — U\ \ — U2\ 

X (-s)mi+m2-f^l -^; -W-l, -"12; S + 1 - mi - 7712; , . (3.2) 

V 1-Vi I-V2J 

To carry out the summations over r and s we employ the integral formula 

Fiia;b,b';c;x,y)= ^f^ /"' ^-i(i _ ^)— i(i _ ^^)-fe(l _ ^y)-^'^^, 

r(a)r(c-a) Jo 

see [m Vol. 1]. In our case b = —r, b' = —s, x = 1 — ^, y = 1 — ^, so in (13. 2p we need to 
compute 

^ <-_-)^ („.„^ (i - , (: - J^)))' (,.. (i - <i - ^)))' (-0.«. 
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^ ( S)mi+m2—i—j — [ ^) ( ^1 '>T'2)mi+m2 

X {1 - (1 - e)(r?ini + r]2Vi) - C{viuiU2 + 772^1 «2)^^+"^~™^""'^ (3.3) 

with the imphcit assumption that rii + n2 > mi + m2- However, by (jl.Sp fa) and (ll.Sp fc) the 
expression in { } vanishes unless ni + 722 = mi + 777-2 • Therefore, ()3.3p becomes 

^ 0mi+m2,ni+n2! 

and consequently, 

Ca = 6m,+m2,n,+n2{l " ^l)'"Hl " ^2)"*^ (r?27^l) 



1, -^^ ) mi +1712 
V^ (-771i)j(-77l2)j / r]lUl Y~^^ / 1 - ttl Y f l-U2 

^ z!j! V 7?2'yiy vi-f^i/ vi-^2 

17 /^ ■ • • , • . W2 . ^^2 

X ti -n2\ -I— 3,1+ J - mi- 7772; -77T.1 - 777-2; 1 , 1 

V Ul Vi 

Smi+m2,ni+n2{l " ^^l)-(l " ^2)^ (7?27;i)-+- ^"^^^ i^l+^^M 



X 



E 



(-777i)j(-7772)j / T^lTil V"'"-' / 1 - Til V/ 1 " 7^2 \ '^ 



z!?! V 772 7; 1 / Vl — 7'l7 Vl — 7i2 



7il7;2 



^ 



X 2F1 -77-2; i+ j - mi- 7772; -7711 - 7T72; 1 , (3.4) 

V ^27^1 / 

by a special case of the last identity (jl.lip . 



4 Final summations in (13.41) 



At the last stage we will set i + j = k, j = k — i, so that the i-sum becomes 

(l-7Xi)(l-7;2)~ 



2-1 



Fi -mi, -k; m2 + I - k; 



{1-U2){l-Vi) 



(-7771 - 7772)fe r^f^ r ^ ^ . -, (1 - 7Xi) (1 - 7;2) 

"27^1 I —mil, —ft, -7771 — 7772, i 



(-77l2)fc V (1 -7i2)(l - 7;i) 

(-7771 -7772)fc „ / , ^ 7Xi7;2\ 

r 2F1 -mi, -k; -mi - 7772; 1 , 

(-"72)fc V 7427^1/ 



by dlS]) and (fTTSjl . Thus 



7-711,712 _ r (la, \'^l(l -,, \'^2(r, -,, \"ii+7?i2 1^771 +7772J. / 772 

'mi,m2 - C'7TM +7712,711+712 1-'- - Vl) [i - V2) (,7?27-'2J 7 TTT 1 

(- jVjmi+m2 V^l 



'^2 



7711+7172 / \"^l/\/;\"2/ I \/ \ i+j 

E [-mi - m2)k sr-^ {-mi)i{-k)i sr^ {-n2,k - mi - m2)j / uiV2\ 
fcl ^—^ iU — m-i — mnV; ^— ' l'(' — m-i — m.n^ ,: \ 



^ /c! ^ ^K-"^i - "i2)i ^ f-{-mi-m2)j V ^^27^1/ 
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since -^ (^l"^) = 1, because r]iui{l - U2) + r?2Vi(l - V2) = 0. But, 

Z^ U {-k)i{k - mi - 1712) j 



now 



fc=o 



k\ 



mi+rn2—i—j 



{-iyi-mi-m2)i+j Yj 

fc=0 

( — 1)*( — "^1 — m2)mi+m2^mi+m2,i+j 
[TTli + m2)mi+m2\^) ^ ^ ^mi+m2,i+j- 



{i + j -mi- m2)k 



k\ 



So 



Tni,n2 
mi,m2 



6m,+m2,n,+n2{l " ^l)™^! " ^2)"^^ ( m^l (l " ^^^ 

{{mi + m2)!)2 ^ (-7ni)i(-l)™i+'"2-^ i 



mi+m2 /^, \ n2 
U2 



U2V1JJ \Ui 

f^2)mi+m2—i 



{-N) 



mi+m.2 —5 i!(-mi - m2)i {-mi - m2)mi+m2-i{mi +m2-i) 



The summand is unless n2 > m,2 + mi — i ^ 712 > m,2, since mi > i. So we set i 
m,i + m2 — n2 + /, / > 0, and get, for the i-sum above 



-"-2)712 (~"^l)mi+m2- 



n2-m2 I ^ 

n2 V^ ("^2 - ■0-2); 



(-mi -m2)mi+m2("T'i +"12)! ^ ^! 



-m2)m2(-"il)mi 



mi\m2l 



{-mi - m2)mi+m2{'mi + m2)! 
since ttt-i + m2 = ni + n2. Thus, 



Jm2,n2 



t I M2 "^2,'^2 

(mi + ?Ti2j! 



7-711,712 

-'mi,r?i2 



57ni,,n2'Jni,n2(l - ^iH (1 - ^2rM - 



^2 



"1 



-r/2i^i 1 



^1^^2 

U2i;i 



X 1 



N 
mi,m2 



mi = ni 



)) 



mi +7712 



To determine the coefficient in terms of rji and r/2, note that 



tilf2 
U2V1 



(1 - Ui)(l - V2) Ui -V1+V2- U2 - {uiV2 - U2V1) 



{1-U2){l-Vi) 

ZJ(l-%-r?2) 

(l_^2)(l-^l)' 



{l-U2){l-Vi) 



from solving ()2.4p . with D = U1V2 — U2V1. 
Now, from ([23]) (a) and (1231) (b) we get 

m = {v2-U2)D~'^, rj^= {ui-vi)D'^, 

while ([r8])(a)-(lLH])(c) give 



m 



Ui Vi 

U1U2 V1V2 
1 — Ui 



1 Vi 
1 flf2 

1 - U2 _ 

fi('(;2 - U2) ^2(^1 - ui) 



1 — ul 



(^2-1) 

Ul(f2-U2) U2{vi-Ui)'' 
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Dr]2Vi{l — vi) 


D 


1 


{l-U2){l-V^) 


V2 — U2 


Vi 


m(l23D. 






Now, 







772^1(1 - V2) = -r/iMi(l - U2), mviil - V2) — = -r/in2(l - U2), 

Ul 

so 

Dr]2Vi{l-V2) U2 _ Di]iU2 _ D _ 1 

{I - U2){1 - Vl) Ul (1-Ul) Vl-Ul %■ 

Thus the normahzation factor is 

(62(7711, m2;iV; %,%)(! -rji -r/2)"^)" • 



5 Proof of the recurrence relation (11.91) 

By the transformation formula (ll.lip . 

P (x. x.)-(l u.rUl ..^x,^ {mi + m2-N)i+,{-m2)k+i{-xi)i+k{-X2)j+i 

-Ul V /^ """1 \'' f ^2 — Ui\ /f2 — 'L'l^ 



1 — til/ \1 — Vl J \ 1 — Ul J \1 — Vl 

So, a straightforward calculation gives 

{N -mi - m2){Pmi+l,m2{xi,X2) - Pmi,m2{xi,X2)) = -(1 - Ui^^ {I - Vi)^^ 

/ d d \ 

^ V^dlZ ^ '"'^d7 ) ^^'^ ~ ^^'^ "^Hl-f^l) ''''Pmum2ixi,X2)), 



where u'^ = Ui/{ui - 1), v'^ = Vi/{vi - 1), i = 1, 2 

ia^ = Ui(l - UiJQi^ 



Clearly u[-^ = Ui{\ — Ui)-Sj-, etc. Hence, with a similar expression for Pmi,m2+i ~ Pmi,m2i 



we can write 

{N — ITli — m2){A{Pm^j^i^rn2 — Pmi,m2) — B{Pm^^rn2+l — Pmi,m2)} 

= B{1 - U2ni - V2r 1^2(1 - U2)^ + t;2(l - V2)^ 

[ OU2 OV2 

X ((1 - 7X2)-"H1 - ^2)-"^Pmi,mJ - ^(1 - Uif^il - ViY^ 

X |ni(l - m)^ + Vi(l - ^l)^} ((1 - ^l)~"Hl - ^l)""^^m,,m.) (5.1) 

for some suitably chosen constants A and B. 

A more convenient form of the right-hand side of (|5.ip is 



{xi{Bu2 - Aui) + X2{BV2 - Avi)}Pmi^m2 + {P ( ( U2 {1 - U2) ^ h f2 (1 - ^2) o" 



9^ + "^^'-"^)9^i 



A i iui{l- Ul) h til (1 - Vl) -— } Pmi,m2- 
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By using the values of A and B indicated in ()1.9p . and those of u's and v's in ()1.3p we find that 

Bu2 - Aui = P1+P2, Bv2- Avi = -{P3+P4). 
Similarly, 



(^Tni[^rmj^ — l^rii2 -'■mi,m2J ^^TT'2{J~^mi,m2 — l ^mi, 



m2) 

What we really need to prove is that 

-(^(1 - ui) + C)ui^ - {A{1 - vi) + C)vi^ 

OUi OVi 

+ {B{1 - U2) + D)U2-^ + (5(1 - V2) + D)V2^ \ Pmum2ixi,X2) = 0. 

Straightforward algebra gives 

A{1 - ui) + C = Pi, A{1 - U2) + C = -P2, 

B{l-u2) + D = -p3, B{l-V2) + D=pi, 

so it amounts to showing that 

d d d d \ 

-PAUI- VP2VI-Z P3^i2 7^ VpiV2-:T— -Pmi/ma = 0. (5.2) 

OUi OVi OU2 OV2J 

Since the w's and v's are expressed in terms of the 4 p's, what we need now is to express the 
derivatives in (15. 2p in terms of those of the p's. Noting that 

{'t' uu -t* Vn -t^ U2 J -^ V2 ) ^ "J X^pu '^P2T ^ps^ -^Pi) > 

for any differentiable function F, with the Jacobian J given by the 4x4 matrix 

/ dui dui dui dui \ 
J = {m„v,„U2„V2,), K.)=^5^,^,9^>5^J > etc., z = l,2, 

we are reduced to the task of proving that 

-PiUiaii + P2Via2i - P3U2a3i + Pi?^204i = 0, (5.3) 

and 3 more similar relations, where the l^l"^^^,- are elements of the inverse matrix J"^, which, 
of course, exists. By a set of long and messy computations we obtain 

an = -piV2A^/p4pIpIS^, a2i = -U2A'^ /pip2plS^, 

asi = -viA^ /pmplS^ a4i = -uiA"^ /p^plS^ (5.4) 



with S = pi + p2 + P3 + Pi, A = piPi — p2P3- Substitution of (15. 4p proves (j5.3p . The three other 
relations are similarly proved. This completes the proof of (jl.9p . 



An additional comment. After this paper was completed we became aware of a recent 
arXiv posting [T3], where the authors point out some important work of H. Mizukawa and 
H. Tanaka [19]. In a future publication we return to the probabilistic origin of the work of 
M. Hoare and M. Rahman and we discuss the relation between the approach in |19j . based on 
the notion of character algebras, and our own. 
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